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4^ The exciting possibility of direct observation of QCD instantons in heavy-ion collisions has re- 

cently been proposed by Kharzeev. The underlying phenomenon, known as the chiral magnetic 
effect, may have been observed recently at RHIC, and a first principles calculation is needed to 

QO confirm and understand the results. The chiral magnetic effect is thought to be visible in the 

symmetric phase, at temperatures above the QCD critical temperature, and in the presence of an 

^— H external magnetic field. We report on first 2+1 flavor, domain wall fermion, QCD+QED dynam- 

^— H ical simulations above the critical temperature, in a fixed topological sector(s), which are used to 

study the electric charge separation produced by the effect. 
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1. Introduction 

Kharzeev proposed that topological charge fluctuations can be observed in the quark-gluon 
plasma created in heavy-ion collisions in the presence of a strong magnetic field [1]. The topologi- 
cal charge fluctuations induce an effective 6 parameter with a spatial gradient which couples to the 
magnetic field to separate positive and negative electric charge[2]. The effect, dubbed the chiral 
magnetic effect, has reportedly been observed by the STAR experiment at RHIC[3, 4]. Because the 
effect depends on the strength of the magnetic field, and its observation on the system being in the 
chirally symmetric phase [5, 6], it is important to study the chiral magnetic effect in the framework 
of lattice QCD in order to confirm this exciting discovery and understand it in detail. The Moscow 
group reported at the meeting on a study in quenched SU (2) gauge theory [7]. 

The physical picture used and referred to throughout this talk is the following. Begin with a 
region where the gluon field strength times its dual, GG, is non-zero. In the classical picture, this 
region corresponds to an instanton, and there will be a chiral zero mode associated with it. This 
zero-mode is an equal mixture of quark and anti-quark components. When the magnetic field B is 
applied over this region, the zero-mode is polarized along the direction of the magnetic field, with 
the quark residing slightly more on one hemisphere and the anti-quark on the other. Thus, there 
is a small polarization of electric charge near the "surface" of the instanton. Because the relevant 
quark degrees of freedom are the low-lying modes of the Dirac operator, it is useful to work with 
its spectral decomposition, and for practical reasons we work with the Hermitian Dirac operator 
B>h = Y5®h> The charge density in units of e in the continuum then becomes 



where the eigenvalues A are real. The extra factor of / is necessary to continue back to Minkowski 
space. A similar expression with similar properties holds in the case of domain wall fermions 
(DWF) at non-zero lattice spacing a[8]. It is easy to see that p = for an exactly chiral mode, so 
that after polarization by the external field, the would-be zero-mode has neither zero eigenvalue nor 
is it exactly chiral. However, it remains nearly chiral as we shall see. Strictly speaking, the lowest 
modes of the DWF Dirac operator are not exact zero-modes, but near zero-modes because of the 
small, explicit, chiral symmetry breaking induced by finite L s . Here we treat this small breaking as 
negligible compared to that induced by the magnetic field, but this still needs to be demonstrated. 

2. Simulation details 

To study the chiral magnetic effect in 2+1 flavor QCD at non-zero temperature and topological 
charge, we use DWF on a lattice of size 16 3 x 8, the RHMC algorithm of Clark and Kennedy, and 
an auxiliary determinant (to fix the topological charge) [9, 10, 11]. The Iwasaki gauge coupling 
is Pqcd = 1-80, extra dimension size L s = 16, light quark mass m/ = 0.013, and strange quark 
mass m s = 0.04. The twisted masses in the auxiliary determinant are £f = 0.0001 and £^ = 0.50 
(see Ref. [9] for definitions). Our simulation parameters were chosen so the system is above the 
critical temperature, T c . To fix the topological charge, the evolution was started with £f = 0.5 and 
gradually reduced to its final value to "freeze" the topological charge to a value of Q = 9 or 10. 
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Though it is a secondary effect, to have the complete electromagnetic picture, we also couple the 
dynamical quarks to photons as well as the external magnetic field. The QED Wilson gauge action 
is used, with coupling Pqed = 1-5. For simplicity for now, the quark charges are set equal. The flux 
of the external field is quantized because of the periodic boundary conditions in units of Injql} 
[12, 13] (q is quark charge), and to maintain gauge invariance, the quarks obey a twisted boundary 
condition [12]. For measurements of the eigenvectors of the Dirac operator, a valence quark value 
of L s = 32 was used in order to obtain low-modes with continuum-like chiral properties (see Fig. 1, 
and for a detailed discussion of the DWF eigenvalue spectrum and how it is computed, Ref. [8]). 




Figure 1: Chirality ( y/j \ Xjfj) of the low-modes of the DWF Dirac operator for the QCD configuration 420 
discussed in the text, with B z = 0. There are 10 zero-modes in this configuration, so the topological charge 
is Q = 10 by the Atiyah-Singer index theorem. 



3. Results 

3.1 Continuum-like instanton 

Before proceeding to the phenomenologically interesting case of 2+1 flavor QCD above T c , 
we investigate the chiral magnetic effect in the background of a single instanton that has been 
discretized onto a lattice[14] with size 8 4 sites. 

-v = -iv"^ (3.D 

p(r) = po ( 1 - — ) S(iWx - r) (3.2) 

\ K max / 

where the continuum instanton of size po/a= 10 is smoothly cut off at a distance from the center 
of r max /a = 3. Because it is cut off, the instanton supports several zero-modes, not just one. In the 
presence of a non-zero magnetic field pointing along the z-direction, layers of positive and negative 
charge appear above and below the instanton in this direction. In Fig. 2, this separation is shown 
for the first zero-mode associated with the instanton. Note, both positive(red) and negative(blue) 
charge appears above and below. However, there is a net positive charge above and minus that 
below so that the total charge on the configuration is zero, as it must be. Whether the dipole 
structure above and below the instanton is a physical effect [15], an artifact of the discretization of 
the continuum instanton, or both, is not clear and needs further investigation. Next, the magnetic 
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flux is doubled, tripled, and quadrupled, and the total charge computed above and below the mid- 
plane of the lattice in the z-direction. The results are plotted in the right panel of Fig. 2. An 
almost linear rise with the strength of the magnetic field is found when only those eigenvectors that 
are nearly chiral, {Xi\T^\Xj) « 1 are used in the charge sums. Higher modes which have chirality 
significantly less than one, but are not paired, are presumably lattice artifacts; paired modes do not 
contribute. The breaking of translation invariance in the x — y plane due to the quark Landau levels, 
which would happen as well in the absence of the instanton[12], is clearly visible in Fig. 2. 



DB: g3-imag,eig.0,eig4.Bx0By0Bz0. 098175. vtk 
Cycle: 98175 
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Figure 2: Left panel: Charge separation computed from a single near-zero-mode for a continuum instanton 
discretized on an 8 4 lattice. B z = 0.098175. Translational invariance is broken in the x — y plane by the 
Landau states of the quarks. Right panel: total amount of charge separated to the lower half of the lattice 
in the z direction for the same configuration. All modes with chirality close to one are included in the total. 
The same amount, but with opposite sign resides in the top half. 

3.2 2+1 flavor QCD 

Having established the chiral magnetic effect in the case of the continuum-like instanton, we 
turn to the relevant case of 2+1 flavor QCD above T c . We study one configuration from our ensem- 
ble, corresponding to monte-carlo time unit 420. Using the 5LI definition of GG [16] and 60 steps 
of APE smearing, the topological charge density for time slice three is shown in the top-left panel 
of Fig. 3. Clearly there are several "lumps" of topological charge and associated localized zero- 
modes of the Hermitian DWF Dirac operator. Note, all of the zero modes (10) are right-handed 
on this configuration and were calculated for a valence quark mass of m v = 0.0001 in order to 
expose their continuum like chiral nature. It is interesting to note that the first non-zero mode pair 
of eigenmodes is localized around an instanton-anti-instanton pair, as shown in the top-left and 
bottom-right panels of Fig. 3. Similar results hold for B z ^ 0. 

Next, we compute the charge density from these zero-modes in the presence of a magnetic field 
in the z-direction, qB z = 2Kn/L x L y , with n = 1, 2, 3, and 4 (q = 1). In Fig. 4 the charge density 
is shown for several values of B z and two different (would-be) zero-modes. Note, even for B z = 
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DB: top.vtk 




Figure 3: Top left panel: Topological charge density. Top right: right-handed zero-mode, eigenvector 0. 
Bottom: right-handed zero-mode, eigenvector 7 (right panel) and eigenvector 1 1 (left). The eigenmodes are 
all localized around the "instantons". B z = 0. 

there is some charge density because the zero-modes are not exactly chiral. It is relative to this 
"background" that we seek the charge separation induced by the magnetic field. The peak values 
appear to be localized around the lumps, possibly suggesting charge is separated above and below 
the centers of topological charge density shown in Fig. 3 (see the bottom panels in Fig. 4). To get at 
the net charge separated on this configuration, the center of the topological charge distribution in the 
z-direction is determined, z = Y^zGG/Y^GG. On configuration 420, z ~ 9.5. The separated charge 
is just the total charge in the two halves of the lattice above and "below" this plane. Unfortunately, 
the electric charge has been computed using the time component of the local current instead of 
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the point-split conserved one. And the total charge on the lattice in the case where B z ^ is not 
conserved compared to the case where B z = 0. This is in marked contrast to the classical instanton 
case described above. Thus we can not conclude that electric charge has been separated. 



DB: q3-imag.eig.420.eigO.BxOByOBzO.vtk DB: g3-imag.eig.420.eigO.BxOByOBzO.0490874.vtk 

Cycle: Cycle: 490874 





Figure 4: Charge density from (near) zero-modes. Top left panel: eigenvector 0, B z — (left) and B z = 
0.0490874 (right); bottom: eigenvector 7, B z = 0.0490874 (left) and B z = 0.0736311 (right). 



4. Summary 

We have reported on a preliminary investigation of the chiral magnetic effect in 2+1 flavor 
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QCD+QED. The effect was clearly established in the case of a classical instanton discretized on a 
four dimensional lattice. On a more interesting 2+1 flavor QCD+QED configuration above T c the 
charge separation around the instantons in an external magnetic field, while suggestive (see Fig. 4), 
was not be clearly observed and needs further investigation. The non-conserved local current was 
used to compute the charge, and the small symmetry breaking induced by finite L s was ignored. The 
results are being updated to address these issues, so we can confirm and understand this exciting 
discovery made at RHIC. In particular, it is important to determine the dependence on temperature 
and the strength of the external magnetic field. It will also be interesting to take full advantage of 
these unique QCD+QED configurations to investigate the correlations of the QED and QCD fields, 
in particular {E • B)q ED (E • B)qcd which we leave for future work. 
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